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We calculate the production of a massive vector boson by quantum fluctuations during inflation. 

This gives a novel dark-matter production mechanism quite distinct from misalignment or ther¬ 
mal production. While scalars and tensors are typically produced with a nearly scale-invariant 
spectrum, surprisingly the vector is produced with a power spectrum peaked at intermediate wave¬ 
lengths. Thus dangerous, long-wavelength, isocurvature perturbations are suppressed. Further, 
at long wavelengths the vector inherits the usual adiabatic, nearly scale-invariant perturbations of 
the inflaton, allowing it to be a good dark matter candidate. The final abundance can be calcu¬ 
lated precisely from the mass and the Hubble scale of inflation, Hi. Saturating the dark matter 
abundance we find a prediction for the mass m « 10 -5 eVx (10 14 GeV/Hi) 4 . High-scale inflation, 
potentially observable in the CMB, motivates an exciting mass range for recently proposed direct 
detection experiments for hidden photon dark matter. Such experiments may be able to reconstruct 
the distinctive, peaked power spectrum, verifying that the dark matter was produced by quantum 
fluctuations during inflation and providing a direct measurement of the scale of inflation. Thus a 
detection would not only be the discovery of dark matter, it would also provide an unexpected probe 
of inflation itself. 
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I. INTRODUCTION 

Inflation mm is a compelling model of the earliest moments of the universe. It addresses many theoretical 
puzzles such as the observed large scale homogeneity and isotropy of the universe. Simultaneously, through 
quantum mechanical perturbations, it seeds density inhomogeneities that can explain the origin of structure 
in the universe [3]. These inhonrogeneities are imprinted on the cosmic microwave background (CMB), and 
where probed, their spectrum matches with those of the predictions of inflation. This remarkable agreement 
is often regarded as the best observational evidence for the inflationary paradigm. 

Cosmological measurements show that the observed growth of structure in our universe requires the 
existence of a new particle, namely, the dark matter. Conventionally, it is assumed that the origins of 
dark matter are decoupled from the mechanics of inflation ■ candidates such as Weakly Interacting Massive 
Particles (WIMPs) are assumed to arise from thermal processes that occur after inflation reheats the universe 
[4] while ultra-light scalars such as axions acquire a cosmological abundance as a virtue of initial conditions 
01S]- While it is reasonable that these two sectors are decoupled, it is tempting to ask if inflation, a theory 
that so beautifully explains the origins of structure in the universe, could also simultaneously be the source 
of the dark matter that is essential for the growth of this structure. 

In fact, inflation contains a natural mechanism to generate a cosmological abundance of particles. The 
transition from the the early inflationary era to the later radiation- or matter-dominated era populates 
field modes that were initially in the vacuum state mm- For example, inflation is expected to produce a 
nearly scale invariant spectrum of gravitational waves [3], whose discovery would be regarded as proof of the 
inflationary paradigm. Similarly, a scalar field whose mass is less than the Hubble scale during inflation will 
generically be coherently populated [3] with a scale invariant spectrum. This could produce an energy density 
equal to that of cosmic dark matter for a sufficiently light scalar (such as an axion). However, the power 
spectrum of such a field would contain isocurvature perturbations that do not match CMB observations. 
This rules out this production mechanism for these dark matter candidates, instead leading to constraints 
on such scenarios in simple models [9]. 

Given these known results for scalars and tensors, we consider the inflationary production of ultra-light 
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vector bosons i.e. vector bosons with a small 100 GeV), but non-zero mass. Such vector bosons may 
arise naturally in frameworks of physics beyond the standard model nni- Further, much like axions, their 
interactions with the standard model can be naturally very small, preventing their thermalization with the 
standard model plasma in the early universe mm- Consequently, an abundance of such particles produced 
during inflation can be cosmologically interesting. We show that ultra-light vector bosons are produced by 
inflation. While they are initially produced with a spectrum similar to that of scalars and tensors, importantly 
their spectrum evolves differently during cosmological expansion. As a result, inflationary production of these 
vector bosons can account for all of the observed dark matter density. Further, unlike the case for scalars, 
we show that the spectrum of density inhomogeneities produced by this mechanism matches with those 
observed in the CMB. This can be accomplished without any additional terms in the Lagrangian (such as 
couplings to the standard model) - all that is necessary is the existence of a massive vector boson. While 
the spectrum is automatically consistent with the CMB, to match the observed dark matter abundance, the 
mass of the vector boson is a function of the inflationary scale. Thus the existence of an ultra-light massive 
vector boson and inflation is sufficient to produce the observed properties of dark matter in the universe. 

Naively, this difference in the spectrum of scalars (or tensors) and vectors would seem to be a surprise. After 
all, the perturbations responsible for the production of these particles is entirely the result of gravitational 
dynamics during inflation. It might seem that the equivalence principle would lead to identical spectra for 
all these bosons. Indeed, for modes that are inside the horizon, the energy density in scalars and vectors 
behaves identically. However, the evolution of the energy density in super horizon modes of the vector boson 
is different from those of scalars. As a result, the final spectrum is different. This is thus consistent with 
the equivalence principle, since the equivalence principle, being a local statement, only guarantees identical 
evolution of sub-horizon fluctuations. The differences in the evolution of super horizon modes is precisely 
the reason for the well known fact that while inflation copiously produces massless scalars and tensors 
(gravitational waves), it does not source massless vectors. 

The cosmological evolution of such ultra-light massive vector bosons have been considered before mm, 
and the results of [12] on the evolution of transverse modes and of m on the evolution of super-horizon 
longitudinal modes agree with our calculations. However, [12; was focussed on preserving a dark matter 
abundance of these vector bosons produced by the misalignment mechanism, wherein the field is generically 
assumed to have a non-zero initial value prior to inflation. Owing to the red-shift in energy density, |12] 
concluded that preservation of the dark matter density through an initial misalignment would require the 
inclusion of additional terms in the Lagrangian with a fine tuned choice of coefficients. Further, they did 
not focus on the fact that a massive vector boson has three degrees of freedom - two transverse modes and 
a longitudinal mode. While the production of the transverse modes is suppressed during inflation (they 
are approximately conformal, much like massless vectors), we find that the longitudinal mode is copiously 
produced. This production is central to our result. The production of the longitudinal mode and its super¬ 
horizon evolution was previously considered in m- However, that paper was focussed on reheating the 
universe through the direct decays of a massive vector boson, and did not consider the late time abundance. 
Since in that scenario the vector was directly responsible for the CMB power spectrum, it needed to be scale 
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invariant, and hence m found that such scenarios are constrained. 

The focus of our work is different - the longitudinal modes of a massive vector boson sourced by inflation 
becomes the dark matter of the universe. Over the wavelengths of perturbations observed in the CMB, the 
dominant density perturbations of the vector field are inherited from the inflaton in the same manner as 
they are for the radiation field, thus matching observations of the CMB. For this to be true, unlike nana, 
we do not need any additional terms - as long as there is a vector boson with the correct mass, inflationary 
production generates a cosmological abundance for it that is consistent with observation. These results m 
demonstrate a novel, purely gravitational production mechanism for dark-matter, which occurs automatically 
during inflation, and appears to be unique to light vector fields. 1 

We begin in section [II] with an executive overview of the mechanism and results. In section III A we derive 
the main result, beginning with the theory of a vector in the expanding universe, following its production and 
subsequent evolution, and finally determining the final relic abundance. The correct imprinting of adiabatic 
density fluctuations is shown in section [TV] We discuss the phenomenology and direct detection prospects of 
such dark matter in section |V| Finally we conclude in section [Vi] 


II. EXECUTIVE SUMMARY 

We calculate the production of a massive vector boson by quantum fluctuations during inflation. The 
longitudinal modes are produced with a non-scale-invariant power spectrum, while the transverse modes are 
suppressed. If the vector is cosmologically stable, the abundance produced ends up as cold relic matter in the 
late universe, in the form of a coherent oscillating condensate of the field. This novel production mechanism 
has interesting consequences that we explore. 

Spectrum For commonly considered light fields, the spectrum of fluctuations produced by inflation is flat 
(or nearly flat) over a large range of wavelengths. This is true for the inflaton itself, for the graviton, and 
for any canonically coupled light scalar field. However, we find that the longitudinal modes of a canonically 
coupled massive vector 2 are produced with a peaked spectrum (see Fig. [4|. The peak occurs at an inter¬ 
mediate wavelength, much smaller than the size of the observable universe but much larger than the usual 
short distance cutoff of the inflationary spectra. The power is greatly suppressed at both large and small 
wavelengths. 

The fact that the spectrum is not scale-invariant arises from an interesting difference in the behavior of a 
scalar field and a vector field in the early universe. The energy density in non-relativistic modes of a scalar 
field is frozen when Hubble is greater than its mass. In contrast, the energy in a vector redshifts as 1/a 2 
in this regime. This affects the long-wavelength modes the most (see Fig. [2]), suppressing power on large 
scales. One important consequence of this is that the misalignment mechanism is ineffective at producing 
a late-time abundance of a vector. Unlike for a scalar such as the QCD axion, energy density stored in 


1 Although for a related production mechanism for extremely heavy dark matter see [IS- 

2 We assume the mass of the vector is “on” during and after inflation, such as in the case of a Stueckelberg (i.e. fundamental) 
mass. Our results do not apply, for example, if the mass is only generated after reheating by the Higgs mechanism. 
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a homogeneous vector field damps away while H > m, making any relic abundance produced in this way 
negligible 3,4 . The difference from a scalar field arises from the fact that the metric (and hence scale factor) 
is inherently part of the norm of the vector (simply relating A M to A M brings in the scale factor). 

Successful density fluctuations This production mechanism is perhaps most interesting as a way to 
generate the measured dark matter abundance. With regard to this, a second important consequence of 
the long-wavelength suppression is the absence of isocurvature modes on cosmological scales, which would 
be visible in the CMB (and are ruled out by observations). For a scalar (such an axion), such isocurvature 
modes are dangerous and typically rule out large parts of parameter space. However, the vector naturally 
avoids these observational constraints. 

Of course, to be a good dark matter candidate the vector must also have the standard scale-invariant, 
adiabatic, Sp/p ~ 1CU 3 4 5 * fluctuations on cosmological scales, which are understood to be the imprint of 
earlier inflaton fluctuations. Luckily, these fluctuations are automatically imprinted onto the density of the 
vector. This can be understood intuitively from the “Separate Universes” approach. Since the observed 
adiabatic fluctuations occur at wavelengths far larger than the scale of the dominant vector fluctuations, 
these wavelengths are far outside the horizon while all the relevant vector dynamics is occurring. Each 
patch of the universe can therefore be treated as a separate homogeneous universe, with the local inflaton 
fluctuation affecting it only as an overall offset in the local clock or scale factor. Because this is just a change 
to the overall clock, it affects massive vector dark matter in the same way as it would affect any type of 
dark matter (e.g. a thermally produced relic particle). So, just as in a normal scenario, the dark matter 
density will fluctuate along with the densities of every other component of the universe. This guarantees the 
correct adiabatic fluctuation spectrum. Thus inflationary fluctuations of a light vector provide a successful, 
and novel, dark matter production mechanism. 

Relic abundance In order to make up the measured dark matter abundance we find a condition on the 
mass of the vector and the Hubble scale of inflation, given by 

n A - ^ cdm x \JQ X 10-6 e y (io 14 Gev) ' ^ 

The current bound on the scale of inflation, Hj <, 10 14 GeV, therefore places a lower limit on the vector 
mass required for this mechanism to generate the entire dark matter abundance, m j> 1CU 5 eV. A lighter 
vector could nonetheless make up an interesting sub-component of the dark matter abundance. Requiring 
m < H during inflation allows the mechanism to work in principle for masses as large as ~ 10 8 GeV. 

Direct detection phenomenology A new vector field will generically have, at some level, a kinetic mixing 
with the photon. This would cause the vector field to couple weakly to charged Standard Model particles, 


3 This applies to both transverse and longitudinal modes, since there is no distinction between them in the zero-momentum 
limit. 

4 The failure of the misalignment mechanism was noted in Ref. m There, a large coupling of the vector to the curvature 7Z 

was introduced to restore the possibility of misalignment production. However, such a coupling would generically be expected 

to introduce a quadratic divergence for the vector mass, and so we do not consider it here. 
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enabling detection of the dark matter with various proposed or existing experimental setups operating in 
different mass ranges [II [Mil. In particular, a relatively high scale of inflation would put the vector in 
the optimal mass range for the experiment recently proposed in m- This experiment would take advantage 
of the fact that the vector field is coherently oscillating at a fixed frequency (set by its mass), to cover a 
wide range of possible vector masses with high sensitivity (see Fig. [6]). 

An interesting feature of the inflationary production mechanism is the large power in fluctuations near 
the peak of the spectrum, at a comoving scale 


-^comoving ~ 10 kill X 


10 -5 eV 


( 2 ) 


On these scales, dark matter would therefore be expected to clump and form self-bound substructures (these 
would form early, resulting in a smaller physical size than indicated in Eq. [2]) . Since the earth is moving 
at around 10 10 km/year, direct detection experiments would be able to map out much of this interesting 
structure on experimental timescales, assuming the structure survives gravitational disruption in the galaxy. 
By both measuring the vector’s mass and inferring its original power spectrum would allow a dramatic 
confirmation of the inflationary origin of the dark matter, as well as giving an entirely new and powerful 
way to probe inflation itself. 


III. RELIC ABUNDANCE OF A VECTOR FROM INFLATIONARY FLUCTUATIONS 

In this section, we show that quantum fluctuations during inflation can produce a calculable relic abundance 
of massive vectors that depends only on the Hubble scale Hj during inflation and the mass m of the 
massive vector boson. The computation is similar to the calculation of inflationary production of scalars and 
tensors. The massive vector field is decomposed into spatial Fourier modes, which are initially sub-horizon 
and evolve to become super-horizon modes during inflation. This evolution populates these Fourier modes. 
Subsequently, the modes re-enter the horizon during radiation domination. To compute the abundance 
and spectrum of the particle today, it is necessary to track the evolution of the mode on both super and 
sub-horizon scales. 

The key point in this calculation is that a massive vector boson has three physical degrees of freedom - 
two transverse modes and one longitudinal mode. The transverse modes, being approximately conformally 
invariant, are not efficiently produced by inflation. The longitudinal mode suffers no such suppression and is 
efficiently produced. This is one of the central results of this section. To show this, we decompose a massive 


vector into transverse and longitudinal modes in sub-section III A and identify the behavior of these modes 
under various regimes. In particular, we identify a regime where the longitudinal modes are relativistic 
and evolve like massless scalars. Using this result, in sub-section |III B[ we show that much like massless 
scalars, inflation will also produce these longitudinal modes. In sub-section |III C[ we track the evolution of 
the produced modes as they transform from being super-horizon to sub-horizon modes. Using the results of 
sub-section III A[ we show that the energy density in non-relativistic, Hubble-damped, super-horizon modes 
of this field redshifts with the expansion of the universe. This is in contrast to massless scalars and results 
in the production of a non-scale-invariant spectrum for the massive vector boson. Finally, in sub-section 
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IIIIDI we calculate the relic abundance of the massive vector as a function of its mass and the Hubble scale 


during inflation. Using this spectrum as the input, we show in section IV that the adiabatic perturbations 
of the inflaton are also imprinted on the spectrum of the massive vector boson, enabling it to be a good dark 
matter candidate. 


A. A massive vector in the expanding universe 

We consider the following action for a massive vector 


s = j V\g\d 3 xdt( - Vie A - l -m 2 g^A,A„) . (3) 

where F is the field strength of the massive vector boson A^. Here, the mass is assumed to be a Stueckelberg 
mass. For an Abelian gauge boson, such a mass term is simply a parameter in the Lagrangian and is 
technically natural. Our analysis will also apply to cases where such a mass term is generated dynamically 
at a scale higher than the inflationary scale. All that is relevant for this analysis is that ([3]) is the correct 
effective action for the massive vector field at all times during and after inflation. This field has three degrees 
of freedom - two transverse modes and one longitudinal mode. Note that the longitudinal mode still exists 
in the spectrum in the limit m —> 0. In this limit, it is decoupled from other degrees of freedom but its 
coupling to gravity is unsuppressed. This is evident from the Goldstone boson equivalence theorem since 
this limit can be obtained by taking the v —» oo, g —> 0 limit of the Higgs mechanism. The coupling of 
the longitudinal mode to gravity is unsuppressed and as we will see below, it is not conformally invariant. 
Hence, it can be efficiently sourced by inflation, much like scalar or tensor fields. 

To analyze these effects, we start with the equations of motion of the massive vector boson field in the 
background of an expanding FRW universe. For the production mechanism, it is sufficient to focus on 
homogenous and isotropic cosmologies as these will generate an initial spectrum for the field. The effects of 
inflationary density fluctuations on this spectrum are phenomenologically important and are calculated in 
section |TV] The metric of an FRW expanding universe is 

ft “' = ( 1 a(t)H J ■ <4) 

where a(t) is the scale factor. In components, the action is then 

S = [ a 3 d 3 x dt-(-^-\d t A — V A t \ 2 —^-|V x A\ 2 + m 2 A? —^-m 2 |A| 2 ^ (5) 

J 2 \or 1 1 a 41 1 a z 1 1 J 

where A = A* (as opposed to than A 1 ). Note that A t does not have a kinetic term and it is thus a non- 
dynamical auxiliary field. Fourier transforming allows us to complete the square and decouple it from the 





physical propagating modes: 


decoupled auxiliary field 


s = 


a 3 d 3 kdt 

(2tt) 3 


1 (k 2 

2fe + “' 


A, - 


ik ■ dtA 


k 2 + a 2 m 2 


— (\ 

2a 2 V 


9t.A\ - 


I k • dtA\ 


- -\k x A\ — rrr \A\ 

n Z I I II 


( 6 ) 


physical modes 


where k is the co-moving momentum. Note that since A(x, t ) is real, its spatial Fourier transform obeys 
A(k,t ) = A(—k,t)*. We can separate A into transverse and longitudinal components At and Al, where 
k ■ A = kAi, and k ■ At = 0. 5 The action separates into an action for the transverse modes, 


^Trails — J 

and one for the longitudinal modes, 

o [ 

^Long — / 


a 3 d 3 kdt 1 


1 (2tt) 3 

2a 2 

a 3 d 3 k dt 

1 ( 

(2tt)3 

2a 2 V 


^ d t A T -^ 


(^2+ m2 )l^| 2 ) - 


d t A L -■ 


At 


O' 


(7) 


( 8 ) 


Switching to the conformal time coordinate dt = adrj , the action for the transverse modes can be written 


St™ = j^\ (Mt| 2 - (<= 2 + aV) |1 T | 2 ) 


(9) 


Here the almost complete disappearance of a displays the well-known conformal symmetry of the transverse 
modes. It is broken only by the small mass term. This means that, unlike the longitudinal modes, production 
of transverse modes by inflationary fluctuations is negligible. We focus on the longitudinal modes for the 
remainder of the paper. 

The action for the longitudinal modes, Eq. (|8j) , has an unusual form. However, field redefinitions allow 
us to recover more familiar forms in the both the high-energy and late-time limits. At early times, when 
am <C k ( i.e. the relativistic limit, where the physical momentum of the mode is much bigger than its mass), 
the a 2 m 2 term can be dropped from the denominator in Eq. (J8|. With a simple field redefinition , 

—* TTf —* 

Tr(k,t) =—A L (k,t), (10) 


the action becomes familiar: 

^Long am<<fc > J df^(|d t 7r| 2 - j^M 2 ) = / a 3 d 3 zdf* ((d t 7r) 2 - ^|Vtt| 2 ) . (11) 

We see that the field i r has the usual action of a real massless scalar field in the expanding universe. As 
expected from the Goldstone Boson Equivalence Theorem, in the highly relativistic limit the longitudinal 
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Note that completing the square to give Eq. [fi] is equivalent to setting At = ik ■ dtA/(k 2 + o?m?) oc Al in the classical 
equations of motion. Assuming At = 0 therefore amounts to throwing out the longitudinal hidden-photon modes. 
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FIG. 1: Cosmological evolution of length scales, and nature of the longitudinal mode in different regimes. The line 
labelled “horizon” shows the comoving horizon size 1/aH , which shrinks during inflation and grows after reheating. 
The line labelled “Compton wavelength” shows the comoving Compton wavelength of the vector, 1/am. Modes of 
the vector field maintain fixed comoving wavevector k , evolving along straight lines from left to right. In the pale 
blue shaded region, where modes are relativistic (m <C k/a), the longitudinal mode behaves identically to a massless 
Nambu-Goldstone boson. In the pale brown shaded region, where Hubble damping is not important (H -C m), 
the longitudinal mode behaves identically to a free massive scalar. In the red triangle between these regions, the 
longitudinal mode has a new behaviour unlike any scalar. Modes crossing the tip of this region reenter the horizon 
just as they become non-relativistic - their wavevector defines the special scale fc*. 


mode is equivalent to a massless (pseudo)scalar. As m —> 0, the interactions of tt with gravity are unaffected 
while it decouples from other sectors of the theory. Hence, this mode can be produced gravitationally even 
when to —> 0. 

At late times, when the Hubble parameter H is small, the expansion of the universe has only a small effect 
on the field. It is then useful to make the field redefinition which makes the longitudinal mode’s kinetic term 
canonical: 


<KM) = 


-A L (k,t). 


\/k 2 + a 2 m 2 

The a-dependence of this redefinition introduces a term of order H 2 into the potential: 

k 2 ' 


, f a 3 d 3 k ,l/i„ , a 2 m 2 rr ,,2 / 2 fc 2 V,| 2 \ 

Long y J2^ d 2\} dt(/)+ k 2 + a 2 m 2H ^ ~ ( W 

j2 dt \ (l 3 ^! 2 “ ( m2 + ^2 ) H 2 + 0 { H 2 )\( t >\ 2 ') 


( 12 ) 

(13) 

(14) 


When to 2 < H 2 , the Hubble term can be dropped (in fact, keeping track of the exact coefficient shows that 
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it can also be dropped at very early times when k 2 a 2 mH 1 but we have already dealt with this regime). 
The action then reduces to the usual action for a real massive scalar in the expanding universe 


c H<g.m x 

*^Long ' 


a^drx dt- (dt4>) -y|V0| - m 


2 i2 


(15) 


In figure[l]we show the different regimes for the behaviour of the longitudinal modes through the expansion 
of the universe. The pale blue region shows the relativistic regime, where longitudinal modes behave the 
same way as modes of a free massless scalar. The late-time behavior is shown in the brown region where the 
modes behave as the modes of a free massive scalar. In this regime, the longitudinal modes are sub-horizon 
and as we will see, the energy density in these modes red-shifts as a -3 . This is of course identical to the red- 
shift of the energy density of matter and this behavior is expected from the equivalence principle. In the red 
shaded triangle, neither of these simplifications apply, and the longitudinal modes have a uniquely vector-like 
behavior which must be determined from its equations of motion. In this regime, the energy density in these 
modes red-shifts differently than the corresponding modes of scalars. However, this behavior does not violate 
the equivalence principle as these longitudinal modes are super-horizon and the equivalence principle is only 
a statement about local interactions of fields in a gravitational background. 

In a universe that undergoes inflation and subsequently reheats, all of these regimes of behavior are 
important. During the inflationary phase, relativistic, sub-horizon modes get stretched, become horizon size 
and eventually exit the horizon. As we will see in sub section fill B[ this process leads to particle production, 
populating the mode. After the end of inflation, the universe reheats and the Hubble scale decreases. These 
modes re-enter the horizon, subsequently evolving to become sub-horizon, non-relativistic modes, during 
which they will go through the red-shaded region. As we will see, the red-shifting of the energy density 
of these longitudinal modes in the red-shaded region is crucial in suppressing dangerous, long-wavelength 
isocurvature power in the spectrum of the massive vector boson. 

The energy density in the longitudinal modes is found using the canonical stress-energy tensor of the field 
which yields p = 2g 00 — £. Fourier transforming the action of Eq. |8jback into coordinate space, and 

restoring g 00 , allows the lagrangian to be written in the form 


2 2 
am 


C 2a 2 ( 9 °° dtAL a 2 m 2 - V 2 
from which the energy density follows, 


d t A L - m 2 A 2 L ) , 


P = 


2 a 2 \ 


dtAr 


2 2 
am 


,d t A L + m 2 A 2 L ) . 


(16) 


(17) 


' a 2 m 2 — V 2 

The Laplacian in the denominator should be understood in terms of its action on Fourier modes, and 
this non-local form arises because we integrated out the time-component of the vector field. In a region 
much longer than the typical wavelength of the hidden-photon field (this applies to all cosmological scales 
of interest), it is appropriate to replace the above expression with its expectation value in that region. We 
can then Fourier transform back into fc-space, and write the energy density as 


M = / 


d In k 


1 


2 2 
am 


2 a 2 \k 2 + i 


d t A L (k, t) + m 2 V AL (fc, t) 


(18) 
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where the power spectrum of a (homogeneously and isotropically distributed) field is defined by 

_ _ _ 2t r 2 

(■ X(k,t)*x(k',t )) = ( 2TT) 3 6 3 (k-k')—Vx(k,t ), 


k 3 


(19) 


so that ( X 2 ) = JdlnkVx(k,t) . 


B. Generation of inflationary fluctuations 

In the previous section, we saw that during inflation, until the modes were well outside the horizon, the 
action for the longitudinal modes is identical to that of a massless scalar field (see figure [I]), up to the 


simple rescaling Eq. (10). This means we can directly apply the standard results for a massless scalar to the 


longitudinal modes of the vector to understand the behavior of these modes during this phase. Light scalar- 
fields are coherently produced during inflation. Beginning their lives as small-scale vacuum fluctuations, 
Fourier modes grow in amplitude as they are expanded beyond the inflationary horizon, and thereafter are 
locked in as classical field fluctuations. 6 . For a free, real, massless, canonical scalar field 7r, it is a standard 
result that modes exiting the horizon are described by 

ik 


n(k,t) = 7r 0 (fc)( 


) ik 

e^JWl 


( 20 ) 


where 7To(fc) is gaussian-distributed with power spectrum (as defined in Eq. (19)) 


Vn 0 (k) = 


Hi 


( 21 ) 


Here Hi is the value of the inflationary Hubble parameter (assumed to be slowly varying) around the time 
the modes exit the horizon. 


Using this result for the longitudinal modes of the vector and using the rescaling Eq. (10), we see that the 
amplitude in the longitudinal mode is 

ik 


^L.exit (M) = A)(£)(l 


a(t)Hj 
2 


e a W H i 


, fkH , V 


( 22 ) 

(23) 


From Eq. (18), we can see that the modes of a given k , at the point when they exit the horizon (i.e. when 
k = aHj), contribute to the energy density as 


dp 

din A: 


(2ny- 


(24) 


6 More technically, the growth of the scale factor a(t) introduces an explicit time dependence into the field’s action. The 
early-time vacuum therefore does not evolve into the late-time vacuum, but into a highly populated state, which can be found 
by performing a Bogolyubov transformation on the initial state. 
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C. Evolution after inflationary production 


In this section, we compute the evolution of the fluctuations after they exit the horizon during inflation as 
well as their evolution from super-horizon modes to sub-horizon modes after the end of inflation. Once the 
modes exit the horizon during inflation, they evolve as coherent classical field modes. The power spectrum 
of the field then evolves as 




A L (k,t) 
A 0 (k ) 


(25) 


where A^ik^t) is the solution of the classical equations of motion with initial condition given by Eq. (22). 
The equation of motion for the longitudinal modes follows trivially from their action: 


( 3, 2 + 


3 k 2 + a 2 m 2 k 2 

k 2 + a 2 m 2 * a 2 


At = 0 . 


(26) 


We now solve this analytically in different regimes. The first three regimes will correspond to the pale blue 
shaded region of figure [T] and we will indeed see identical behavior to a free massless scalar. The fourth 
regime we consider corresponds to the brown shaded region of figure [l] and there as expected we will see 
matter-like behavior. Finally we turn to the red shaded, “uniquely vector-like” region of figure [l] where we 
will discover new behavior unlike that of a scalar. We summarize the results in figure [2] 

• Horizon exit during inflation H = Hi to. We can drop to 2 , and the equation of motion is 
approximately 


(d 2 + 3 Hjd t + k 2 / a 2 ) A l » 0 


( ad a ad a + 3ad a + k 2 /{aHi) 2 )A L « 0 . 


(27) 


The solution to this a term of the form Eq. (22), plus second term whose form is the complex conjugate 
of the first. We see that the inflationary initial condition (i.e. de Sitter vacuum) gives rise to only the 
first term, setting the second to zero. 

• Super-horizon relativistic regime - H k/a to. This is the region labelled “frozen kinetic 
energy” in figure [ 2 ] Here we can drop the k 2 /a 2 and m 2 terms, and the equation of motion is approx¬ 
imately 


(d 2 + 3Hd t )A L » 0 


(d a aHd a + 3 Hd a )A L « 0. 


(28) 


During inflation, when H = Hi ss const, this has solution Al = ci + C 2 a -3 , while during radiation 
domination, when H oc a~ 2 , this has solution Al = c\ + C 2 <i~ 1 . In either case the second term dies 
fast, and so in this regime 


At, = A n = const 


(29) 


to an extremely high accuracy. The contribution to the energy density from these modes therefore 
evolves as 


2 a2 / 2 

p ~ m A L /a oc a 


( 30 ) 
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FIG. 2: Evolution of the energy density in longitudinal modes, from inflationary production through to matter 
radiation equality. As expected from figure [l] the evolution is the same as it would be for a massive scalar in all 
regions except the red triangle labelled “Vector regime”. In that regime, the energy stored in a scalar would be 
constant, whereas for the vector it damps as a -2 . This damping suppresses large-scale isocurvature modes, allowing 
the produced vector abundance to make up the dark matter. This abundance is dominated by modes of comoving 
size 1 /k t indicated by the dashed line. The details of reheating do not affect these modes, as long at it occurs before 
they reenter the horizon. 


This behaviour is the same as it would be for a massless scalar, whose kinetic energy damps as a 2 in 
this regime while the field itself is frozen. This is just as expected from the discussion in section [III A| 

• Sub-horizon relativistic regime - k/a rn. //. This is the region labelled “radiation” in hgure|2] 
Here we can drop the m 2 terms, and treat the Hd t term as a small perturbation. Switching to conformal 
time dt = adi], the equation of motion is approximately 

( d 2 + 3Hd t + k 2 /cl 2 )Al ~ 0 >—k (<9“ + k~ + ‘2aHd r ^')A x ~ 0, (31) 

with solution 

A l ~ 1 (c 1 e ikr > + C2 e- ikr ‘). 
a v y 


( 32 ) 
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The energy density in these modes evolving as radiation, 

p ~ m 2 A 2 L /a 2 oc a~ 4 , (33) 

which again is the same as the well known behaviour of a massless scalar in this regime. 

• Late time non-relativistic regime - m^$> k/a, H. This is the region labelled “matter” in figure [2] 
Here we can drop the k/a terms, and treat the Hdt term as a small perturbation. The equation of 
motion is approximately 

(<9 2 + Hd t +m 2 )A L « 0, (34) 

with solution 

A L ^ ^( Cl e imt + c 2 e- imt ) . (35) 

V a 

The energy density in these modes evolves as matter 

p ~ m 2 A|/a 2 oc a -3 , (36) 

just as it would for a massive scalar in this regime. 

• Hubble-damped, non-relativistic regime - H 3> m k/a. This is the region labelled “vector 
regime” in figure [ 2 } We can drop the k 2 /a 2 and to 2 terms, and the equation of motion is approximately 

(d 2 + Hd t )A L « 0 ^ (d a aHd a + Hd a )A L « 0. (37) 

During inflation, when H = Hj « const, this is solved by A^ = C\ + C 2 <x _1 . As in the super¬ 
horizon relativistic regime, we can drop the rapidly-decaying second term. However, during radiation 
domination, when H oc a~ 2 , the solution is 

A l = ci + c 2 a (38) 

Now there is a growing term, which one might guess would dominate the solution. However, this regime 
follows a long period of the super-horizon relativistic regime, in which the field became constant to 
an extremely high degree. The continuity of A and d a A then prevents the linearly-growing term from 
taking off in the current regime, and the correct solution is again 

A L = A 0 = constant. (39) 

This result is not entirely obvious, and needs checking with care - we do this in appendix [B] The 
evolution of the energy density in these modes damps as a -2 . This is quite different to the energy 
density of massive scalar modes in this regime. In that case, the field <j> would be frozen, contributing 
to a constant vacuum energy m 2 4> 2 . 

! m 2 A 2 L /a 2 a a~ 2 vector 
m 2 4> 2 — const 


scalar 


( 40 ) 
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FIG. 3: Primordial power spectrum of the amplitude of a massive vector field’s longitudinal modes, produced purely 
by inflationary fluctuations. The spectrum is shown at a time just 3 e-folds after H = m, and shorter wavelength 
modes (on the right of the plot) are still relativistic. At later times the fc -1 scaling will continue all the way to the 
right of the plot. 


We can now see that the power spectrum of the longitudinal modes of vectors produced by inflation is 
different from those of scalars. For both the vector and the scalar, modes of all values of k have the same 
energy density when they exit the horizon, dp/dink ~ Hj. The subsequent redshift of this energy density 
varies for different modes and at different times as indicated in figure [2] In particular, the energy density in 
vectors and scalars redslrifts identically for modes that are in the pale blue and brown regions of figure [2] 
This is unsurprising since in these regions, the effective action of the vector is identical to that of scalars 
(with some field redefinitions). However, in the red-region of figure [2j as seen in (40), the energy density in 
vectors redshifts as a~ 2 while the energy density in a scalar would be constant. 

For the vector, it can be seen fairly simply (geometrically) from figure [2] that the modes whose energy 
redshifts the least are those of the special wavenumber fc*. Modes with longer wavelengths (further up the 
figure) receive extra redshifting in the “vector regime”, causing the final energy in them to fall as k 2 at small k. 
Modes with shorter wavelength (further down the figure) undergo rapid redshifting while they are behaving 
like radiation, causing the final energy density in them to fall as k at large k. This means that the power 
spectrum of the vector has a peaked structure (shown in figures [3] and [4]) , with power concentrated at the 
wavenumber /c*. This implies that the power in experimentally probed (CMB) long wavelength fluctuations 
of the energy density of vectors is suppressed. This production mechanism does imply significant power at 
short wavelengths, but these have not yet been accessed by experiment. 

This is in stark contrast to scalars. The energy density in a scalar mode is constant when it is in the 
red-region of figure [2] Hence, all modes with momentum lower than fc* have the same power, leading to 
isocurvature power at length scales probed by the CMB. 
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D. Final relic abundance 


We can now make use of the scalings described in the previous subsection, and calculate the final abundance 
of cold vector matter that was produced. As shown in figure [3j the power spectrum of the massive vector 
is peaked at the wavenumber fc*. This means we can estimate the final abundance by considering only the 
contribution from k ss A;*. The energy density starts as Hf/(2ir) 2 at horizon exit, redshifts as a ~ 2 until 
a = a*, and then redshifts like matter, a* is the value of a when H = m, and = a*m. These are 
approximately related to Hubble at matter-radiation equality by 


tt* — d\H—r 


//„ 


x a n 


fc* = a*m « y/Ha 


,m x a 


m.r.e. • 


This gives an estimate of the abundance at matter-radiation equality 


Hf 


Pve ctor 


or, using H 1 


2 

m.r.e. 


Pcdm(d 

n 


( 2^) 2 

)/M\ 


/ ^exit \ ^ / d* \ ^ 

V a* / V a / 

/Mpn 


I 2 PI 


HfH&.r.e.mi 


(27T) 2 


^ a ln . r . e . ^ ' 


T 




Hjiriz 


f^cdm 


(2t r) 2 M|,T m . r . e . 

To get a more precise result, we can combine Eqs. (18, [25) to write the vector abundance as 


VP? f / Q? \ 

Pvector(A) = ^ 2 J ^ ( ^2 _|_ a 2 m 2 '^ >d t A l (^) *) + T’Al {k, t) J 

fc 2 


Hj a*fc 2 


( a 

A L (k,t) 

2 ^ 
a 6 

d t A L (k,t ) 

V 

Ao(k ) 

a*(/c 2 + a 2 m 2 ) 

Ao(k) 


(41) 

(42) 

(43) 

(44) 

(45) 


87t 2 a 3 J 

Here t ) is the solution to the classical equations of motion, with initial conditions at horizon exit given 


by Eq. (22). The integral is an 0(1) factor which can be found by numerically solving for Ap(k,t). At 
late times it is independent of m, and we find it to be approximately equal to 1.2. The yjm scaling comes 
from the a*fc 2 factor, which can be calculated precisely from the definitions above. Given the observed dark 
matter abundance, p c dm(to) = 1-26 x 10 -6 GeV/cm 3 , we find the final result 

ftvector _ / / Hj \ (HPi 

O cdm ~ V 6 X 10-6 eV l 10 14 GeV J ' 1 ’ 


E. Comment on misalignment production 

Previous work has considered production of a vector dark-matter abundance from cosmological initial con¬ 
ditions nu da 1 19) . Our results are different from the results of this previous work. The reasons for the 
differences are as follows. References mm added a specific coupling to the scalar curvature in the la- 
grangian, while we do not include this term. Reference m showed that, in contrast to the case of a scalar, 
misalignment is ineffective at producing a vector abundance, unless this large coupling to the scalar curva¬ 
ture is added. We agree with this result. Further we note that in this case, inflationary fluctuations cannot 
produce the dark matter, so one must use misalignment. There is in general both misalignment production 
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and production by inflationary fluctuations. When the coupling to the scalar curvature is added so that the 
massive vector behaves as a scalar, inflationary fluctuations will produce the normal, flat spectrum, just as 
for a scalar. However if this makes up all of the dark matter it is ruled out by the bound on isocurvature 
perturbations by many orders of magnitude. Thus in the case where the coupling to the scalar curvature is 
added, misalignment production (or some other mechanism) must dominate over production from the infla¬ 
tionary fluctuations by orders of magnitude. Additionally, we note that the large coupling to the curvature 
is not consistent with the treatment of m 2 as a small spurion of gauge symmetry breaking, and would be 
expected to introduce a quadratic divergence to the mass. 

As noted in na, without the large coupling to the scalar curvature, misalignment production is ineffective 
and cannot generate the dark matter abundance. This also applies to any long wavelength mode produced 
by inflationary fluctuations, as can be seen from the low k part of Figure [3j We suspect however that it 
may be possible to have something like misalignment production of a massive vector, in a different way. 
The norm of the vector decreases during inflation as we have shown, ex a -2 . Thus if we want a 

uniform (misalignment produced) massive vector field to be dark matter today, then at some point during 
inflation this norm would have passed very large scales (e.g. the Planck scale), and the energy density also 
would have been large (larger than the total inflationary energy density). This is another way to state why 
misalignment production does not work in the standard case. However it is reasonable to suspect that the 
effective field theory of this massive vector breaks down at some point. So beyond some field excursion 
the potential may change from a simple quadratic term. In fact it may not be described as a simple, 4D 
field theory at all - for example it may arise from a higher dimensional string theory. Such a degree of 
freedom may naturally redshift very differently than a massive vector during inflation and we would need 
the full UV completion. Clearly then the initial condition of when during inflation the massive vector starts 
from its effective cutoff is unknown but would determine the final abundance and may allow misalignment 
production. This is a change to the massive vector at some point in its field space. There may also be a 
dynamical change to it in the early universe either during inflation or radiation dominance. If the massive 
vector is produced from some other field (e.g. the way axions can be produced from strings) then this may 
create a dark matter abundance of the vector in a very different way which is not subject to these objections 
to misalignment production. Thus, while misalignment production does not work for the canonical massive 
vector field considered here, we note for completeness that there may well be other production mechanisms 
which make other predictions. 


IV. ADIABATIC AND ISOCURVATURE DENSITY FLUCTUATIONS 


In order to be a good dark matter candidate, the massive vector must have the correct, observed power 
spectrum of the matter density of our universe. In particular, it must not have large isocurvature perturba¬ 
tions on long length scales. This is normally an issue for a light, bosonic field, but as we discuss in Section 


IV A the massive vector naturally avoids this problem. Further, it must have the observed adiabatic, nearly 


scale-invariant fluctuations on cosmological length scales. We demonstrate in Section [IVB] that this is indeed 
the case. 





18 



FIG. 4: (Lower plot) Primordial density power spectrum of a massive vector produced by inflationary fluctuations. 
The spectrum is shown at a time when all modes covered are non-relativistic, but before self-gravitation of the modes 
is important (this should be approximately valid until matter-radiation equality). The peaked part of the spectrum at 
large k is the isocurvature power produced by the inflationary fluctuations of the field itself. The also flat part of the 
spectrum at small k corresponds to the usual adiabatic fluctuations, which are imprinted onto the field from inflaton 
fluctuations. The values m = 10 -5 eV and Hi ~ 10 14 GeV were used here, corresponding to fc* « 1400 pc^ 1 . Because 
the density fluctuations are 0(1) for k ~ fc*, fluctuations on shorter scales (in the gray region) are not well described 
by the power spectrum alone (see Fig.[5|. These higher k modes are also expected to be affected by quantum pressure 
later in their evolution. For comparison, the top plot shows the power spectrum of the field amplitude (figure [3]). 


A. Isocurvature Fluctuations 

As demonstrated in Section [III| the power spectrum of the massive vector field is peaked as in figure |3j with 
0(1) power at the special scale /c*. Since they are generation without correlated fluctuations in the radiation 
bath, these are isocurvature fluctuations. Isocurvature perturbations are dangerous, since on cosmological 
scales are strongly constrained by CMB observations. However, this limit only applies to scales within a few 
orders of magnitude of the current size of the universe. On the other hand, K corresponds to a cosmologically 
tiny scale, 

1 /K ~ 10 10 km x v/ 10 ~ 5eV . (47) 

V m 

The isocurvature power spectrum fall off on long wavelengths (low k), to unobservable levels. To demonstrate 
explicitly this we calculate the density power spectrum in appendix [A] The result is plotted as the peaked 
spectrum in the lower panel of Figure [4] The power falls as k 3 below K, becoming completely negligible on 
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FIG. 5: Typical variation of the field and energy density along a line. This is a random example, generated according 
to the power spectra shown in Fig. |4j i.e. before structure formation occurs. (Note the y- and ^-components of the 
field are not shown, and so the lower plot is not exactly the square of the top plot.) As can be seen, the density is 
dominated by “lumps” of size L ~ 7r /k*. On longer scales, overdensities are caused by random clustering of these 
lumps. On smaller scales, there is complicated sub-structure both within and between the main lumps. If the vector 
is discovered in direct detection experiments, its present-day field profile will be mapped out as the experiments sweep 
though the dark matter halo. 



cosmological scales. We conclude that isocurvature perturbations are not a problem for the massive vector 
produced by inflation. 

To get a better picture of what this field looks like, one random realization of short-wavelength fluctuations 
of the vector field with the corresponding energy density is shown in Figure [5] The 0( 1) fluctuations of the 
field and energy density are clear. This could have interesting consequences for formation of dark matter 
structures on these length scales, as well as for direct detection experiments. This is discussed further in 
Section [V| 

B. Adiabatic Fluctuations 

Intuitively the massive vector will pick up the adiabatic fluctuations of the inflaton because all of the 
abundance is generated by small wavelength fluctuations (small compared to cosmological scales). Thus 
on much longer scales it looks homogeneous. This is just like a WIMP, QCD axion, or any other dark 
matter candidate. The distribution of WIMPs appears homogeneous on long scales but in fact there is large 
inhomogeneity at scales around the average inter-WIMP spacing. Of course this does not affect the fact 
that the WIMP picks up the adiabatic fluctuations of the inflaton. Although the peak in the massive vector 
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power spectrum (found in Section III) is at much longer scales than typical inter-WIMP spacings, it is still 
very short compared to cosmological scales. So we know that the massive vector will indeed pick up the 
adiabatic fluctuations of the inflaton and so can be a good dark matter candidate. In this subsection we 
show that this intuitive argument is correct. 

We wish to calculate the spectrum of density perturbations imprinted on the massive vector field by the 
fluctuations of the inflaton. The relevant inflaton fluctuations that are observed in the CMB are much longer 
wavelength than the modes of the massive vector that carry the dominant power. This large separation of 
scales allows us to use the “separate universes” approximation. We will calculate the energy density of the 
massive vector assuming a homogeneous inflaton perturbation to the metric across the entire universe. We 
follow the conventions of Dodelson (20 and work in Newtonian Gauge and, as usual, ignore anisotropic 
stress. The perturbations are then described by substituting g 00 —► g 00 (l + 2 $) and a 2 —► (1 + 2 <f>)a 2 in 
the unperturbed metric. The perturbations $(fc,t) are zero during inflation, then turn on around reheating, 
and remain constant until they re-enter the horizon in the late universe. Since $ varies over wavelengths 
given by the size of the fluctuations seen in the CMB, and since these wavelengths are so much longer than 
the massive vector modes, $ can be taken to be homogeneous. Then we compare the energy density in 
different “universes” (actually just different large regions of the universe) with different values of $. In this 
sign convention adiabatic perturbations to the matter density scale as p oc (1 + 3 $). By deriving this answer 
for the massive vector we will show that the fluctuations of the massive vector on the long length scales 
observed in the CMB are indeed adiabatic, with the correct, nearly scale-invariant spectrum arising from 
inflation. This allows the massive vector to be a good dark matter candidate. 

Consider a region (or “separate universe”) with a metric perturbation. We can find the action for the mas¬ 
sive vector in this background by using the standard action for the longitudinal modes Eqn. ([8]), reinstating 
g 00 with the perturbed value —(1 + 2 $), and making the replacement a 2 —> (1 + 2 $)a 2 . Thus, working only 
to linear order in <f>, we find 

r a 3 (l + 3 $)d 3 kdt 1 — 2$/ (1 + 4 $>)g 00 a 2 m 2 


L ° ns J (27r) 3 \/|<7 00 1 (1 + $) 2 ° 2 V/c 2 + (1 + 2$)a 2 m 5 

This can be simplified by defining new variables 

k' = k( 1 - 2$) 
m' = m( 1 — $). 


I r\ \ I 2 21/1 I 2 

\d t A L \ - m 2 \A l \ 


(48) 


(49) 

(50) 


Then the action becomes 


^Long — 


a 3 (l + 3$)d 3 fc dt 1 / 
(27t) 3 ^H(1 + $) 2a 2 V 


2 /2 
a z m 


k' 2 + a 2 m' 2 


d t A L 


AM 2 )- 


(51) 


Note the similarity with the usual action for the unperturbed massive vector Eqn. (with $ = 0). In fact, 
the equations of motion are exactly the same as for the unperturbed massive vector except with k and m 
replaced by k! and nn !, since the overall constant out front is not relevant for this. 

Thus we know how the massive vector will evolve in this “universe.” It will create a peaked power spectrum 
as shown above. However the wavelength of the peak will change slightly, as will a*, the scale factor at which 
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the massive vector started acting like matter (as in Fig. [2j. We now determine a* by the scale factor when 
H = m!. Since H oc a~ 2 during radiation dominance, we can relate this to the scale factor of matter-radiation 
equality a mre 


a* 

^mre 


h, 


= ^ = (i + $) 

m 


//„ 


m 


(52) 


where H mie is the Hubble constant at matter-radiation equality. Note that matter-radiation equality is still 
defined globally to be the same time for all the “separate universes” (so we do not put a prime on a mre ). This 
is because we want to know the energy density at a single globally defined time across the entire universe. 
We are only concerned with the $ dependence of our answers. And the peak of the power spectrum for the 
massive vector occurs at k' = fc* which is defined by the crossing-point in Figure jlj so we have = ^ 7 . 
This gives us a peak at 

k * = flmreV kI mre Ul f = (1 2^)^ mre ^mre^- (53) 


Of course we want the actual wavelength at which the peak occurs, in other words we want the value of k 


not the value of k' at which the peak occurs. These are related by Eqn. (49), so the peak is at 

3 

k = k* = (1 -}- 24>)/c^ = (1 T ^ T) o Tnre \JH mre m. 


(54) 


These formulae determine the power spectrum of the vector in this “separate universe” with nonzero 4). 

We also need to know the amplitude of the massive vector field produced during inflation in this universe. 


The power produced during inflation is given by Eqn. (23), but now we need to modify it to take into account 


the nonzero <f>. We will only be interested in the power at the peak so we will use k given in Eqn. (54). This 
does change with 4> because the wavelength of the peak in this “universe” is physically different. However 


in Eqn. (23) we must use the unperturbed mass m instead of m!. This is because, in Newtonian Gauge, 4> 
is zero up until reheating where it the jumps to its constant value for the rest of the history of the universe. 
Thus 4> = 0 during inflation and so the power produced around the peak wavelength is 


^Ao(M = 


k,H, y _ 
27 im J 


n 2 TAT tr2 


(55) 


(2ir) 2 m 

After inflation the value of the field A in the longitudinal mode with k = k* is then fixed until the crossing 
point a*, as in Figure]!] Note that one could worry that this statement is no longer true at reheating because 
$ jumps rapidly there. However one can show from the equations of motion that this rapid change in $ will 


not change the value of the field. Thus from Eqn. (551 we have the power in the dominant field mode when 
it begins to act like matter. 

We are interested in the total energy density as a function of $. This is derived in the same way as 


Eqn. (18) was derived for the unperturbed universe. But note that we want the energy density as seen by a 
global observer, not an observer in each “separate universe.” Thus we use the global definition of position 
(and wavevector) so we Fourier transform in k not k' to find the energy density. In the action for the 


longitudinal modes, Eqn. (51), the first factor in the integral is just the usual integral over phase space (in 
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global coordinates). Fourier transforming that action back into position space then gives an energy density 


m = / 


d in k 


1 


2 a 2 V k' 2 + ( 


j'Pdt a l (k, t ) + m' 2 V AL ( k , t ) 


(56) 


which is the same as Eqn. (18) except the fc’s and m’s inside the integral have been changed to k' and m! 
(though note that the argument of the power V is not changed). In the regime where the vector acts as cold 
dark matter, the kinetic energy and the potential energy are the same on average by the virial theorem. So 
we can just take twice the second term for the average energy density: 

1 


P(t) = J c 


dinfc — (1 - 2$)m z VA L {k,t) 


(57) 


Finally the power at matter-radiation equality can be found in terms of the power produced by inflation 

once it is acting as 




Eqn. (55) by redshifting the power by a factor since the field redshfits oc 
matter (as in Fig. [2]). We will focus only on the power in the dominant mode so the integral is removed. 
Using Eqns. (l52| and (f55l) in Eqn. (57) then gives the energy density at matter-radiation equality 


i + 5® 




777-2 




(58) 


This is the correct result for adiabatic perturbations, that the energy density scales as (l + §$). All of our 
sign conventions can be checked by following the above procedure for a scalar instead of a vector, in which 
case one finds the same factor (l + §<f>). Thus the massive vector dark matter does indeed pick up the nearly 
scale-invariant, adiabatic perturbations of the inflation. This is illustrated in the lower panel of Figure [4] So 
the massive vector can be a good dark matter candidate. 


V. DARK MATTER PHENOMENOLOGY 

We have shown that inflation automatically generates a cosmological abundance for any massive vector 
boson, with the abundance completely determined by the Hubble scale during inflation and the mass of the 
vector boson. For masses ~ 6 x 10 -6 eV (lO 14 GeV/Hj) 4 , this abundance is equal to the observed dark 
matter density. For lower masses, it becomes a subdominant component of dark matter. Given the model 
independent nature of this mechanism, there is a strong case to search for this abundance through laboratory 
experiments. Detection in the laboratory requires the massive vector boson to interact non-gravitationally 
with the standard model. 

It is reasonable to expect such interactions. For example, the massive vector boson could kinetically mix 
with the photon, though the operator 

C D \eF^F^ , (59) 

where here F’ llv represents the new vector’s field strength. As a dimension four operator allowed by symmetry, 
it is unsuppressed by heavy mass scales and is generically expected to exist, although different UV completions 
allow e to be naturally extremely small ED- Such an interaction can lead to a wide variety of possible 
constraints and new searches - see for example pza ns na Eaazi and other references given below in this 
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FIG. 6: Prospects for direct detection of a cosmological vector abundance through kinetic mixing with the photon. 
Direct detection experiments are sensitive to the combination £y r py shown on the y-axis (my is the mass of the new 
vector, and py its cosmic abundance - this plot does not assume the vector makes up all the dark matter). In the 
vertical pink band, high-scale inflation (10 13 GeV <J Hi <, 10 14 GeV) generates the full dark matter abundance for 
the vector. Masses to the right of this band require lower values of Hi, while to the left of this band the vector can 
constitute an interesting dark matter subcomponent. The blue dashed line shows the projected reach of a recently 
proposed search with resonant LC oscillators na. with the blue dotted line showing the potential improvement by 
multiplexing in the high-frequency regime. Tan shaded regions are excluded due to cosmological disturbances caused 
by the relic vectors, while green regions are excluded by null results from the XenonlO and ADMX dark-matter 
searches. Gray-blue shaded regions indicate constraints on the existence of the vector independent of its cosmological 
abundance (these are truly constraints on e, so to plot these here it was assumed that the vector makes up the largest 
possible fraction of the dark matter that is consistent with inflationary production). See section |v] for more details. 


section. If so, over a wide range of masses 10 12 — 10 3 eV, a cosmological abundance of such particles can 
be detected using well developed electromagnetic resonator technologies as suggested in [HCEH1 HE]. The 
estimated reach of such experiments is plotted in figure [6] These experiments are sensitive to the combination 
£ \J~Py where py is the local energy density in the massive vector boson, which of course, cannot be larger 
than the dark matter density. Since the sensitivity scales as yjpy, and since e is a priori unknown, these 
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detectors have the ability to probe vector boson densities that are a small subcomponent of the dark matter 
density. This is particularly important in this scenario since the abundance generated by inflation will be 
smaller than the dark matter density for vectors with mass less than ~ 10” 5 eV (lO 14 GeV/fi/) 4 . Given our 
ignorance of the ultraviolet physics responsible for the origins of these massive vector bosons and the model 
independent way in which the vector boson abundance is generated, there is a strong case to search for such 
a cosmic abundance over a wide range of parameters. 

Interestingly, as is clear from figure [6j electromagnetic resonator technology is particularly well suited 
to search for massive vector bosons in a several-decade mass range around ~ 10 -6 eV, corresponding to 
resonator frequencies ~ 10 8 Hz. If the scale of inflation is close to the present experimental bound of 
Hi ~ 10 14 GeV, a vector boson in this mass range will have an abundance equal to the observed dark matter 
density. There is thus strong motivation to leverage the existence of high precision electromagnetic sensors 
in this frequency range as a positive detection could yield strong evidence for inflation. The blue dashed 
line shows the projected reach of a recently proposed search with resonant LC oscillators |T5], with the blue 
dotted line showing the potential improvement by multiplexing in the high-frequency regime. 

As is evident from figure [6j there are a variety of constraints on such kinetically mixed vector bosons. 
A cosmic abundance of these vector bosons is constrained in the tan shaded regions. The central tan 
shaded region in figure [6] is ruled out by CMB distortions due to conversion of the vector into photons (T2j. 
The rightmost region by is excluded due to depletion of the vector by conversion into photons. We derive 
this bound following the analysis of m, but taking the initial abundance to be the maximum allowed 
by inflationary production, whereas in [12] an exponentially large initial abundance was considered. The 
leftmost, semi-transparent tan region is claimed in [12] to be excluded due to depletion of the vector by 
conversion into plasma modes. This analysis appears to have ignored the (potentially large) disruptive effect 
this would have on the baryon plasma, and so may not be valid. We do not attempt to reanalyze this 
constraint here. As pointed out in [19] . conventional dark matter direct detection experiments are sensitive 
to vector dark matter if its mass is above the experimental threshold energy. Similarly, as pointed out in m, 
the ADMX axion dark matter search is also sensitive to vector dark matter. The green shaded regions in [6] 
show the corresponding exclusions from results from the ADMX [39] and the XenonlO m experiments. 

The blue-gray shaded regions are excluded regardless of the cosmic abundance of the vector boson - they 
only rely on the existence of the vector boson as a degree of freedom in the theory. These bounds include 
stellar production of the vector mm, precision tests of electromagnetism [29] [30] [53], and distortion of the 
CMB due to conversion of photons into the vector [53] . Since these are bounds on e rather than Ey/py , an 
assumption about the vector abundance was needed in order to put the bounds on the plot. In the vertical 
pink band and to its right, the full vector dark matter abundance can be generated by inflation consistently 
with bounds on the inflationary scale {Hj <[ 10 14 GeV), and we assume this to be the case. To the left of 
this band, we use the maximum abundance that can be generated by inflationary fluctuations (p oc y/m). 
With this assumption, the black region in the upper left would require e > 1 which is forbidden. 

The primordial dark matter power spectrum produced by this mechanism (see figure [4]) is peaked at 
comoving scales of order k~ 3 ~ 10 lo km x y/l0~ 5 eV/m, and drops as k 3 for larger comoving scales. Thus, 
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over a range of length scales near this peak, the density perturbations are significantly bigger than the 
adiabatic density fluctuations of ~ 10~ 5 . They can therefore be expected to go non-linear and become self 
bound before large-scale-structure (e.g. galaxy) formation, resulting in significant small-scale dark-matter 
substructure. At the peak of the power spectrum, the density fluctuations are 0(1), and hence we expect 
these fluctuations to become self-gravitating around matter-radiation equality. These peak fluctuations would 
then decouple from the subsequent expansion of the universe, leading to structures in the dark matter density 
today at distances ~0.1AU, with over-densities of ~10 5 GeV/cm -3 . Longer wavelength fluctuations, with 
lower power, can be expected to decouple somewhat later, leading to a range of dark-matter over-densities 
today on length scales not too far from 1AU - i.e. roughly the distance the earth moves in a year. Structure 
on these scales is therefore traversed by laboratory experiments over experimentally accessible time scales. 
Remarkably, direct detection experiments should therefore be able to directly probe this structure, measuring 
a range around the peak of the power spectrum. This would give dramatic evidence for this production 
mechanism, enabling a direct probe of inflationary dynamics through the physics of dark matter. For these 
ideas to bear fruition, it is necessary to understand the non-linear formation and effects of galactic dynamics 
on such structures. It would also be interesting to see if these structures have other observable astrophysical 
effects. We will explore these consequences in future work. 

VI. CONCLUSIONS 

We have shown that the observed dark matter abundance and the spectrum of its density inhomogeneities 
can naturally be produced by inflation if the dark matter particle was an ultra-light vector boson. The 
production mechanism is entirely gravitational it only requires the existence of a massive vector boson 
without significant self interactions. As long as the coupling between the standard model plasma and this 
vector boson is small, the abundance of the massive vector boson produced by inflation will not thermalize 
and will contribute to the dark matter density today. Unlike scalars, inflationary production of massive 
vector boson dark matter does not lead to large scale isocurvature perturbations on the CMB. Instead, the 
dark matter density has large isocurvature perturbations only at short distances that have not been probed 
by experiment. On the large distances probed by the CMB, it has suppressed isocurvature power and its 
power spectrum is dominated by the usual adiabatic density perturbations that are imprinted on it from 
inflaton fluctuations. 

The abundance produced by this mechanism is calculable - it is completely determined by the Hubble scale 
during inflation (specifically, when the dominant modes exit the horizon) and the mass of the vector boson. 
This is in contrast to other production mechanisms that are often discussed for ultra-light bosonic dark matter 
candidates such as axions. For example, the density produced by the misalignment production mechanism 
[II 0 entirely depends upon unknown initial conditions. Moreover, the misalignment mechanism does not 


7 We estimate that field-gradient pressure may be significant for modes shorter than k* , possibly preventing their collapse, 
but that longer scale modes will behave as cold pressureless matter, allowing them to form structures. 
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generate a significant cosmic abundance of massive vectors unless there are large additional interactions. 
Unlike scalars, the energy density in the massive vector boson field redshifts rapidly during inflation, diluting 
the initial abundance in the field. While this dilution can be avoided through the introduction of additional 
finely tuned interactions |12j . or if the misalignment was generated near the end of inflation, they add 
complexity to this scenario. Other mechanisms such as the emission of these particles from topological 
defects such as strings may yield calculable abundances but are beset with computational challenges 1551 . 
Further, it may be difficult to obtain direct observational evidence for such events in the early universe. 
If the Hubble scale during inflation is inferred through observations of the primordial gravitational wave 
spectrum, our mechanism predicts the approximate mass (to « 6 x 10” 6 eV (lO 14 GeV/FZ/) 4 ) a vector boson 
must have in order to constitute all of the dark matter of the universe. 

If the inflationary Hubble scale was in the range 10 13 —10 14 GeV that is accessible to the next generation 
of primordial gravitational wave experiments, the predicted range of the vector boson dark matter mass 
would be ~10 -5 -10 -1 eV. These vector bosons may interact with the standard model through kinetic 
mixing, enabling the possibility of direct detection. Of course, in contrast to WIMP dark matter, since the 
inflationary production is gravitational, such interactions are not essential for this mechanism. However, if 
such interactions exist, massive vector boson dark matter can be searched for over this range of masses using 
well developed electromagnetic resonator technologies, as suggested in ,18, 48]. Further, the existence of 
new vector bosons in this frequency range can also be directly probed by various “Light-Shining-Through-a- 
Wall” experiments wherein these vector bosons are produced by a laboratory source and then subsequently 
detected USES]. Thus, in this range of parameters, this scenario leads to an exciting set of experiments. If 
vector boson dark matter is detected in the mass range 10 -5 —10 _1 eV, it would give strong impetus to the 
next generation of primordial gravitational wave detectors. 

This mechanism, being gravitational, will generate a cosmological abundance for a massive vector field, 
even when its mass is not 6 x 10~ 6 eV (lO 14 GeV/FZj) 4 . For lower masses, the abundance generated will be 
subdominant to the dark matter density. It is important to devise experiments to search for such light vector 
boson dark matter even if they are a subdominant component of the dark matter density, since a discovery 
in such a channel provide a new experimental probe of inflation. Conversely, if a high inflationary scale is 
experimentally observed, our results would seem to constrain larger vector masses. However, it should also 
be noted that our results are subject to assumptions about the behavior of the horizon around the time 
the dominant modes exit and reenter it. For example, the universe could have been reheated to a very low 
temperature following a period of matter domination. In that case, the dominant vector modes would begin 
to redshift like matter earlier, resulting in a lower final abundance, and allowing a significantly larger vector 
mass. 

The power spectrum of the dark matter density produced by this mechanism is peaked at comoving scales 
k- 1 ~ lO 10 km x yGO -5 eV/m. Dark matter perturbations at this comoving scale become self-gravitating 
during the era of matter-radiation equality and decouple from the subsequent expansion of the universe, 
leading to structure in the dark matter density at distances ~ 0.1 AU. It is thus possible to sample different 
parts of this structure over ^year time-scales, and hence to directly measure it in the laboratory with direct- 
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detection experiments. The detection of such short distance structure in the dark matter would give dramatic 
evidence for this production mechanism, and a would provide a new probe of inflation itself. 

A central reason for the success of this production mechanism is the fact that the dark matter density 
inherits large distance adiabatic perturbations. Hence, even if the dark matter production mechanism has 
significant power at short distances, as long as the mechanism does not directly produce similar power at 
large distances, it will reproduce the observations of the CMB. It would be interesting to see if there are 
other classes of theories where such a suppression may naturally occur. These might enable the existence of 
structure in the dark matter at length scales ( e.g. ~1 AU) that are short on galactic length scales but large 
enough to be relevant for laboratory experiments. 

Cosmic Inflation is theoretically compelling. Weakly interacting, ultra-light fields such as the graviton are 
naturally produced during inflation and as a consequence of their suppressed interactions, their abundance 
is not thermalized by the plasma. Gravitational waves are known to exist and are hence a natural target 
for experimental searches. However, in light of the challenges that must be overcome to detect them, it is 
interesting to ask if inflationary dynamics can be probed through the existence of other weakly interacting, 
ultra-light fields. Such fields exist naturally in many frameworks of physics beyond the standard model 
and they can interact with the standard model in a variety of ways, unlike the restricted interactions of 
gravitational waves. The remarkable advances in the field of precision metrology over the past two decades 
have made it possible for us to search for a cosmological abundance of these particles over a wide range of 
masses and interactions. A discovery in such an experiment would probe the ultra-high energy scales that 
are responsible for the production of these particles, the mechanics of inflation and the subsequent evolution 
of the universe. 
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Appendix A: Power spectrum of the Held and the density 

The power spectrum Vx{k,t) of a (homogeneously and isotropically distributed) field X is defined here as 

_ _ 2tt 2 

(X(k,t)*X(k',t)) = (2tt f6 3 (k-k')—Vx(k,t), (Al) 

so that 

JdlnkV x {k,t) = (X 2 ). (A2) 

At late times the energy density p{x) in the hidden photon field is proportional to |A(£)| 2 . Defining the 
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density fluctuation S(k) in the usual way, 

p(x) = (p) (l + J -0^S(k)e^ , (A3) 

S(k) is given by 

*>■**>■ (A4) 

If all the hidden photon inodes are longitudinal, satisfying A(k) = kA^ik), then this becomes 

(A5) 

Since the primordial distribution for A^k) is gaussian, we can find the power spectrum of 5 in the following 
way: 
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Integrating over <f> and changing variables from cos 6 to p 
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— q\ gives the simple formula 
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(AS) 


Using the results of section m we calculate this for the vector field at late times. For small k, Vs falls as 
k 3 , while at large k it falls as k _1 . The full numerically calculated power spectrum is shown in figure |4j 
Note that the distribution of the density fluctuations are not gaussian. At long wavelengths, where the 
fluctuations about the mean are small, gaussianity is a good approximation. However at the special scale 
&*, fluctuations are large and the distribution is non gaussian. On shorter length scales, the locally defined 
mean density varies dramatically from the global mean density, as can be see in Fig. [5j On these scales it is 
not obvious that S(k) is even a useful quantity to consider, since it is defined globally. 


Appendix B: Careful treatment of the super-horizon regime 

In section [III C| we concluded that when H is the largest scale in the equations of motion, the solution to 
the field equation is for A^(k, a) to remain constant. This is in spite of the fact that, for k < fc*, once modes 
enter the non-relativistic post-reheating regime, their equations of motion have a solution that grows linearly 
with a until a = a*. The reason this growing term can be dropped is that, prior to entering this regime, the 
equations of motion have already driven the a-derivative of Ai(/c,a) to be extremely small. We now check 
this fact more carefully. 















29 


We know that Eq. (22) is the correct solution before and during horizon exit, up to possible corrections of 


order to/ Hi which could be present because m = 0 was taken when deriving Eq. (221. Shortly after horizon 


exit, Eq. (22) behaves as 


_ _ / k 1 2 ik 3 

AUk,«)-+Mk)( 1 + 2^5 + 3^3 
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(post-horizon exit). (Bl) 


We therefore write An(k, a) = Ao(k)(l + e(fc, a)), with e behaving, shortly after horizon exit but still during 
inflation, as 
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We then solve the equation of motion for e and check that it remains small. 
Using dt = aH(a)d a , the equation of motion for e can be rewritten as 
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e has been dropped on the right hand side, since we are assuming it to be small and checking for consistency. 
We know H « Hj during inflation, and H « m(a*/a) 2 as a approaches a*. We can therefore approximate 
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(It can be checked that deviations from this that occur well after horizon exit and well before a = a* do not 
affect the conclusion.) 


The general solution to Eqs. (|B3| |B4| is 
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The terms in the first line dominate at the beginning of the period, and can be matched onto the desired 


behavior in Eq. (B2), giving C\ = 0(rn/Hi) and C 2 = 1 + 0(m/Hi). The terms in the second line grow and 
dominate at the end of the period, near a = a*. The first of these is the linearly growing term that we were 
concerned about. Having determined that C 2 ~ 1, we now see that this term is indeed incredibly suppressed. 
Note that the remaining two terms, which are not suppressed by to/ Hi, are simply the 0(1) corrections that 
are expected as the solution transitions into the matter-like regime at a = a*. 
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